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ABSTRACT 

The relation between ices in the envelopes and disks surrounding YSOs and those in the quiescent 
interstellar medium is investigated. For a sample of 31 stars behind isolated dense cores, ground-based 
and Spitzer spectra and photometry in the 1-25 /zm wavelength range are combined. The baseline 
for the broad and overlapping ice features is modeled, using calculated spectra of giants, H2O ice 
and silicates. The adopted extinction curve is derived empirically. Its high resolution allows for the 
separation of continuum and feature extinction. The extinction between 13-25 /xm is ^^50% relative 
to that at 2.2 /im. The strengths of the 6.0 and 6.85 /im absorption bands are in line with those 
of YSOs. Thus, their carriers, which, besides H2O and CH3OH, may include NH|, HCOOH, H2CO 
and NH3, are readily formed in the dense core phase, before stars form. The 3.53 ^m C-H stretching 
mode of solid CII3OH was discovered. The CII3OII/II2O abundance ratios of 5-12% are larger than 
upper limits in the Taurus molecular cloud. The initial ice composition, before star formation occurs, 
therefore depends on the environment. Signs of thermal and energetic processing that were found 
toward some YSOs are absent in the ices toward background stars. Finally, the peak optical depth of 
the 9.7 i^m band of silicates relative to the continuum extinction at 2 .2 nm. is significant ly shallower 
than in the diffuse interstellar medium. This extends the results of iChiar et al.l (|2007l ) to a larger 
sample and higher extinctions. 

Subject headings: infrared: ISM — ISM: molecules — ISM: abundances — stars: formation — infrared: 
stars — astrochemistry 



1. INTRODUCTION 

Planets and comets are assembled from the gas, dust 
and ices in circumstellar disks that in turn evolve from 
the envelopes and clouds surrounding YSOs. Many key 
questions about the molecular composition of protostel- 
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lar environments are still unanswered. The infrared spec- 
tra of YSOs are rich in absorption features caused by the 
vibrational transitions of silicates and of solid H2O, CO, 
i^CO, CH3OH, NH3, OCS, CH4, CO2, "CO2, an d pos 
sibly of OCN~, NHt, HCOOH and HCOO ~ (Gibb et al.l 
2OOI iBoogert fc Ehrenfreundl l2004: Boog ert et al.l 1200 



Oberg et al.l 120081: iZasowski et~l.l i2009i: iBottinelli et al.l 
20101 ). Band profile studies revealed that both CO 



and CO2 are present in multiple components of mixed 
or pure ices that are formed or destroyed at dif- 
ferent stages of the cloud and envelope evolution 
iTielens et al."1991':'Gerakines et al."1999 HBoogert et al.l 
J002; Pontoppidan ct al. 2003a, 2008). It was also shown 
that the prominent 6.0 and 6.85 fim absorption fea- 
tures toward low and high mass YSOs are caused by 
the overlapping modes of simple sp ecies: H2O, H2CO , 
HCOOH, NH 3, NH +, and CH3OH (IKeane et al.ll2001bl: 



IBoogert et al.|[2008l ). Not all absorption observed in this 
wavelength region is explained, however, such as the very 
broad component labeled "C5" in Boog ert e t al. (2008). 
Ices processed by heat, UV photons or cosmic ray hits 
are p otential candidates. Both ions (iSchutte fc Khannal 
I2003f) and hydrocarbons (jCibb k Whittetll2002D have ab- 
sorption features that are spread over much of the 5-8 fim 
region. PAH molecules embedded in the ices could be re- 
sponsible for some of the absorption as well. 

These molecules are formed by a complex inter- 
play between numerous formation and destruction pro- 
cesses. Grain surface chemistry likely dominates the 
formation of molecules (e.g., H2O, CH4, CO2) in cold, 
dense molecular clouds, leading to icy mantles (e.g., 
Tielens_& Hagcn 1982). In the protostellar environment 
these ices may sublimate through the influence of shocks 
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or photons in the inner envelopes and disk surfaces. The 
formation of complex molecules in the ices may be trig- 
gered by heating and by t he absorpti on of energ etic pho- 
tons (iSchutte et alJ Il993t iGrccnbcr ^ et all [1995) origi- 
nating in the protostar-disk boundary layer. In the most 
shielded areas (outer envelopes, disk mid-planes) deeply 
penetrating cosmic rays can induce similar effects. It is 
still quite uncertain which molecule formation and de- 
struction processes are important in these different en- 
vironments. What is the composition of ices in pristine 
clouds, prior to star formation? To what degree does 
the solar system molecular inventory resemble quiescent 
cloud inventories? 

A rigorous observational approach is required to fur- 
ther constrain the origin of unidentified absorption 
bands, especially in the 5-8 /im region, and to disentangle 
the many effects of the physical environment on molec- 
ular evolution. An established technique is to contrast 
the molecular composition toward YSOs with that to- 
ward backgrou nd stars tracing quie scent molecular cloud 
material (e.s... iWhittet et al.iTl983l [l998). Large abun- 
dances of the most volatile species (e.g., pure, "apo- 
lar" CO ice) are observed in quiescent lines of sight, 
which is attributed to the absence of nearby heating 
sources. Likewise, amorphous ice structures are observed 
in these sight-lines, whereas toward YSOs they are some- 
times crystalline (Smith e t al.l Il989t iGerakines et al.l 
ll999HPontoppidan et al. 2000?^ Tight upper limits have 
been set to the abundance of the sp ecies OCN~ (cyan ate: 
IWhittet et aDr2001.) and CH3OH (IChiar et al.lll99l to- 
ward background stars, suggesting that specific physical 
conditions found only around YSOs may be needed to 
initiate their formation. 

So far, mid-infrared (A > 5 ^m) studies have focused 
on only a few background stars: four toward the Tau- 
rus molecular cloud, one t oward Serpens, and one to- 
ward IC 5146. iKnez et al.l ()2005D show that the absorp- 
tion bands in the 5-8 /Ltm wavelength region are similar 
to those toward YSOs. The 15 /xm CO2 band shows 
a non-polar ice component consistent with the pristine 
nature of the Taurus cloud (|Bergin et al.l [20051 ). a con- 
clusi qn that can be draw n for Serpens and IC 5146 as 
well (IWhittet et al.l [20091) . 

The studied samples are too small to fully assess the ice 
inventory prior to star formation, given the wide range of 
environmental conditions possible. Much larger samples 
of background stars are needed, covering wider ranges of 
extinctions and cloud types. Selecting these samples is 
now possible using the catalogs produced by wide field 
Spitzer IRAC and MIPS mapping projects, in combi- 
nation with the Two Micron AU Sky Survey (2MASS; 
iSkrutskie et al.l [2006( ) . In this work, background stars 
are selected from the catalogs of nearby (<350 pc), iso- 
lated, compact (< 5 arcmin), dense (> 10** cm~^) cores 
of the Spitzer Legacy team "From Molecular Cores to 
Planet Forming Disks" fc2d: lEvans et al.ll2003[ ). These 
cores were originally identified in optically selected re- 
gions of extinction by observing transi tions of NH3 sen- 
sitive to high densities (> 10"^ cm~^; iMvers fc BensorJ 
|1983() . Based on IRAS colors, a subdivision was made 
between starless cores and cores containing stars, and 
it wa s concluded that rou ghly half the cores contain 
stars (jBeichman et al.lll986D . With the sensitive Spitzer 



MIPS and IRAC cameras, deep searches revealed new 
YSOs, down to stellar masses of 0.01 Mq, even in 
some cores that we re previously categorized as "starless" 
([Young et al.|[200l . 

Studying the interstellar medium in isolated cores is 
attractive, because it lacks the environmental influences 
of outflows and the resulting turbulence often thought 
to dominate in regions of clustered star formation (e.g., 
p Oph, Serpens). Their star formation time scales can 
be significantly larger due to the dominance of magnetic 
fields over turbulenc e and the result ing slow process of 
ambipolar diffusion (|Shu et al.l[l987HEvanslll999( ). The 
longer time scales may alter the gas composition, such 
as reduced C/CO and H/H2 ratios, which would qualita- 
tively alter the grain surface chemistry. Furthermore, the 
longer exposure time of the ices to cosmic rays or cosmic 
ray induced ultraviolet photons may initiate reactions 
with large energ y barriers leading to the formation of 
comp lex species |Greenberg et al.lll995t iBernstein et al.l 
[2OO2I) . 

While ices in the quiescent medium toward isolated 
dense cores have rarely been studied, t he 9.7 /xm band 
of silicates received more attention. Chi ar et al.l (|2007f ) 
found that toward the IC 5146, B68 and B59 cores, 
as well as some lines of sight through the Serpens and 
Taurus molecular clouds, the peak optical depth of the 
9.7 fim band relative to the K-band extinction is signifi- 
cantly shallower compared to the diffuse medium. Here, 
this relation will be revisited with a much larger sample. 
With a simultaneous analysis of the ices bands, some 
possible explanations of this relation can be addressed, 
e.g., the effects of grain growth. 

Also, in the process of disentangling ice and dust fea- 
tures from the attenuated stellar spectra, an extinction 
curve needs to be assumed. Broad-band and spectro- 
scopic studies have sh own that the extinction remains 
high even at 25 /Ltm ([Chapman et al.l [20091 IMcClurd 
I2009D . In this work, this result will be re-assessed, using 
Spitzer spectra and source-specific models. In the pro- 
cess a high resolution extinction curve will be produced, 
in which absorption features and continuum extinction 
can be separated. 

Finally, this study of isolated dense cores is part of a 
larger program to investigate dust and ices in a range 
of environments. While here many cores are observed 
with a few background stars per core, IChiar et al.l ([20TT[ ) 
study many sight-lines behind one specific core in the 
IC 5146 region in much detail. Efforts are under way to 
study large nearby clouds as well. 

The selection of the background stars is described in 
Sj2l and the reduction of the ground-based and Spitzer 
spectra in iJ3| In ij4. li the procedure to fit the stellar con- 
tinua is presented. This is a crucial step in the analysis. 
For this, a "high resolution" extinction curve is needed 
in which ice and silicate features can be separated from 
continuum extinction. This extinction curve is derived 
in M.2\ Subsequently, in the peak and integrated 
optical depths of the ice and dust features are derived, 
as well as column densities for the known identifications. 
Then in ij4.4i the derived parameters are correlated with 
each other and with previously studied YSOs. In H5.ll 
the comparison of ice abundances in the background stars 
and YSOs is discussed, followed by a discussion on ob- 
servational tracers of ice processing in H5.2[ H5.3\ shows 
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how the observations of the ices further constrain expla- 
nations for the deviating extinction curve and Ak versus 
r9.7 relation in dense clouds. Finally, the conclusions are 
summarized and an outlook to future studies is presented 
in gi 

2. SOURCE SELECTION 

Background stars were selected from a sample of 
69 isolated molecular cloud cores that were mapped 
with 6'iPzteer/IRAC and MIPS by the c2d Legacy team 
(jEvans et al.ll2003L l2007f ). The cores are limited by size 
(< 5 arcmin), distance (< 350 pc), and density (> 10** 
cm^-^), and both historically desig nated starless and star- 
forming cores (from IRAS colors; iBeichman et al.lll986D 
are included. The catalogs generated by the c2d team 
were systematically searched for sources meeting the fol- 
lowing criteria: 

1. The overall SED (2MASS 1-2 ^m, IRAC 3-8 ^im, 
MIPS 24 ^m) must be that of a reddened Rayleigh- 
Jeans curve. 

2. Regions of high extinction {Ay > 5 mag), associ- 
ated with the core, must be traced. The regional 
extinction was measured using the N ear-Infrared 
Color Excess (NICE) method (e.g., iLada et all 
[1999), which uses five 2MASS sources nearest each 
target to derive statistically reliable extinctions. 

3. Fluxes must be high enough (>5 mJy at 8.0 //m 
or >30 mJy at 15 /xm) to obtain Spitzer/IRS spec- 
tra of high quality (S/N>100) within ^30 minutes 
of observing time per module. This is needed to 
detect the often weak ice absorption features and 
determine their shapes and peak positions. 

Thus, a sample of 32 candidate background stars was 
selected. One of those (2MASS J17111631-2725144 be- 
hind B59) will not be analyzed here, as its IRS spectrum 
shows a silicate emission band and no ice bands. The 
remaining 31 sources (Table [1]) span a range of A^ val- 
ues and are located behind 16 isolated cores. Eight out 
of sixteen cores are "starless" in the IRAS definition, 
although some of tho se have gjpzteer-det ected YSO can- 
didates (e.g., L1014: lYoung et al.ll200l . The selected 
background stars are usually within ~60" of these YSO 
candidates. Spitzer spectra of many YSO candidates 
were obtained, but they will be presented in a forthcom- 
ing paper. 

Finally, to compare the analy sis methods pre - 
sented here with previous work (jKnez "eTall [2005h . 
two well-studied stars behind the Taurus and Ser- 
pens molecular clouds were included in the sample: 
2MASS J04393886-f 261 1266 (Elia 3-16) and 2MASS 
J18300061+0115201 (CK 2, [EC92] 118). 

3. OBSERVATIONS AND DATA REDUCTION 

Spitzer/mS spectra of background stars toward iso- 
lated dense cores were obtained as part of a dedicated 
Open Time program (PID 20604). The complemen- 
tary Taurus and Serpens background stars were ob- 
served in the c2d program (PID 172). Table [1] lists all 



sources with their AOR keys, and the IRS modules they 
were observed in. The spectra were extracted and cal- 
ibrated from the two-dimensional "BCD" spectra pro- 
duced by the standard Spiteer pipeline (version S15.3.0), 
using exact ly the same method and routines discussed in 
iBoogert et al. (2008). Uncertainties {la) for each spec- 
tral point were calculated using the "func" frames pro- 
vided by the Spitzer pipeline. 

The Spitzer spectra were complemented by ground- 
based Keck/NIRSPEC (jMcLean et all [1998^ H, K and 
L-band spectra at a resolving power of i?=2000. These 
were reduced in a way standard for ground-based long- 
slit spectra, using bright, nearby main sequence stars 
as telluric and photometric standards. In the end, 
all spectra were multiplied along the flux scale in 
order to matc h near-infrared broad-band photometry 
from 2MASS (Skrutski e et all 120061 ) as well as Spitzer 
IRAC and MIPS photometry from the c2d catalogs 
(jEvans et al.ll2007D . using the appropriate filter profiles. 
The same photometry is used in the continuum deter- 
mination discussed in ||4l Catalog flags were taken into 
account, such that the photometry of sources listed as be- 
ing confused within a 2" radius or being located within 
2" of a mosaic edge were treated as upper limits. The c2d 
catalogs do not include flags for saturation. Therefore, 
photometry exceeding the IRAC saturation limits (at the 
appropriate integration times) were flagged as lower lim- 
its. Finally, since the 2MASS and Spitzer photometry 
represent different surveys, and their relative calibra- 
tion is very important for this work, the uncertainties 
in the Spitzer photometry were increased by adding the 
zero-point magnitu de uncertainties listed in Table 21 of 
lEvans et all (|2007l) . 

4. RESULTS 

The observed spectra (left panels of Fig. [T]) show many 
distinct absorption features on top of reddened stel- 
lar continua: 3.0, 6.0, 6.8, 9.7, and 15 ^m. These 
are attributed to ices and silicates, an d were previ- 
ously observed toward niany Y SOs (e.g.. IBoogert et all 
120081 : iPont oppidan et al.' '2 0081) and some background 
stars (e.g., Knez et al. 2005 nBergin et al.ll2005l ). At a 
lower level, weak features from the stellar atmosphere 
are present as well (e.g., 2.4 and 8.0 /im). The separa- 
tion of interstellar and photospheric features is essential 
for this work and is discussed next. 

4.1. Continuum Determination 

The continua and some of the absorption features in 
the background star spectra were fitted by minimizing 
the reduced-x^ of a model consisting of the following pa- 
rameters: 

• Spectral type. An accurate spectral type determi- 
nation is necessary to separate interstellar features 
and features from the stellar atmosphere, as well 
as to know the overall spectral energy distribution. 
Correction of the stellar CO and SiO bands at 5.3 
and 8.0 ^m is particularly important, as they blend 
with the 6.0 /im ice band and the 9.7 /im band of sil- 
icates. One-dimensional, opacity-sampled spectra 
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Table 1 

Source Sample 



Source 


Core'' 


AOR keyt^ 


Module^ 




d 


2MASS J 










-loouzyy 


IHAM 04191 


0014900992 


SL 




1 A 


UOUOZioO— 


-oyuyou4 


BHR 16 


or\i yioncjioo 

0014896128 


SL, LL2 






UoUyoioO— 


-oDU4UoO 


UCjj dU-ol 


00148984^2 


SL, LL2 






UoUyo4Do— 




r^^^ Qn Qi * 


nni A onci 7ft 


GT 






1 oni /I Qni 


-DOUo4zz 


T^f 007 V O fi* 


nm /ion7i Qft 
UUi4yu / ioD 


GT T T 






lzU14oyo— 


-DOUoOoD 


T\f 007 V O Q* 
DK^ Zy i . (-Z.O 


nm /ion7Qoo 
uui4yu / oyz 


GT T T 

Oij, LiLiZ 






1 c;/i oi C/i v 


-OZ4o14d 


T\r^ Q070 1 1 c 

JJU oz / Z-f-lo 


uui4oyy4oD 








io4ziDyy — 


-0Z4 / 4oy 


T^r^ Q070_Ll fi 
oZi Z-f-lo 


uui4oyyyDo 


GT 






1 i iiioUi- 


0*70^; ion 


"D c;n* 


nm ACir\r\AQr\ 
UU14yUU4oU 


GT 


2.41-4 


14 


1 1 lilooo- 


070 '/I /I /I 

-Z( Zi 144 


r> oy 


nm /I Qop; 1 r\A 
UU14oyoiU4 


GT 


9/11 A 


1 A 

14 


1 "71 1 onnc; 
1 t liZUUo- 


07071 Q 1 


13 oy 


nm A QC\A QAQ 

UUi4oy4o4o 


GT 


9 QQ A 

z.oo-4 


1 A 

14 


T "71 c; c; c '70 

1 ^ iooo i o- 


on c; c; Q 1 o 


L 100* 


or\i yinr\i r\ 

0014901760 


OT 


2.38-4 


14 


1 1 idU4d (- 


on 7n7o 


L 100* 


or\i /mr\or\T 

0014902016 


OT 

bL 


2.39-4 


14 


171dU»oU- 


on cr o 1 ,10 

-2Uool42 


L 100* 


nm /inm o/ic 
UUi4yUi24o 


GT 


2.38-4 


14 


1814(J71z- 


r\7r\o /I T o 

0708413 


L 438 


OOI49O0O08 


GT 

bL 


2.82-4 


14 


18160600- 


-0225539 


CB 130-3 


0014896897^ 


SL 


2.82-4 


14 


18165296- 


-1801287 


L 328 


0014908928 


SL 


2.41-4 


14 


1 Qi fic;ni v 
iSlDOyl 1 - 


1 Qm 1 t;Q 
-ioUlioo 


L 328 


0014907648 


SL, LL2S 


9/11 A 

z.41-4 


1 A 

14 


18170426- 


-1802408 


L 328 


0014907904 


SL, LL2S 


2.06-4 


14 


18170429- 


-1802540 


L 328 


0014904064 


SL, LL2 


2.39-4 


09 


18170470- 


-0814495 


L 429-C 


0014908160 


SL 


2.38-4 


14 


18170957- 


-0814136 




nm A fioso/1 /I 
uui4oyoy44 


GT T T 


2.38-4 


14 


18171181- 


-0814012 


L 429-C 


0014908416 


SL 


2.06-4 


14 


18171366- 


-0813188 


L 429-C 


0014908672 


SL 


2.06-4 


14 


18172690- 


-0438406 


L 483* 


0014905600 


SL, LL2 


2.38-4 


14 


192015974 


-1135146 


CB 188* 


0014897408 


SL 


2.41-4 


14 


192016224 


-1136292 


CB 188* 


0014897664 


SL, LL2 


2.06-4 


14 


192144804 


-1121203 


L 673-7 


0014904576 


SL 


2.06-4 


14 


212405174 


-4959100 


L 1014 


0014902528 


SL, LL2 


1.49-4 


14 


212406144 


-4958310 


L 1014 


0014909696 


SL 


1.49-4 


14 


220637734 


-5904520 


L 1165* 


0014903040 


SL 


1.49-4 


14 


043938864 


-2611266'^ 


Taurus MC 


0005637632 


SL, SH 


2.59-5 


3lj 


183000614 


-0115201' 


Serpens MC 


0011828224 


SL, SH 


2.82-4 


14 



^ Star-forming cores as defined by IRAS studies are indicated with a . 
^ Identification number for Spitzer observations 

Spitzer/YBS modules used: SL=Short-Low (5-14 pm, _R ~ 100), LL2=Long-Low 2 
(14-21.3 Aim, R ~ 100), SH=Short-High (10-20 ^l^a, R ~ 600) 

Wavelength coverage of complementary near-infrared ground-based observations, 
excluding the ranges ~1. 79-2. 06 and ~2. 55-2. 82 ^m blocked by the Earth's atmo- 
sphere. 

" Averaged with AOR key 0014899712 

^ Also observed in AOR key 0014896896, but rejected due to high solar activity 
^ LL2 observed, but not included in the analysis because of source confusion 
^ Well-studied background star Elia 3-16 JKnez et al.lf2005l and references therein) 
' Well-studied background star [EC92] 118 (IKnez et al.lf2005l and references therein) 
J Near-infrared data obtained from the ISO archive dWhittet et al.lll998f) 



were calculated usin g the MARCS 2 model atmo- 
spheres following De cin et al.l ()2004l ). Only giants 
(luminosity class III) were considered (Table [2]), as 
these are bright, and most likely to have made it 
into the source selection. The blend of CO overtone 
lines at 2.25-2.60 /im is a sensitive tracer of spectral 
type. For the later spectral types (>M0 III), the 
4.7 /im fundamental CO transitions become promi- 
nent as well. Although most of this region was not 
observed, its edge at 4.0-4.1 /im was; it drops off 
sharply for spectral types >M0 III and was used as 
another independent indicator of spectral type. For 
the earliest spectral types (late G), the 2.16 /im H-I 
Br7 line is important. Finally, because the opacity- 
sampled model spectral lines can only be compared 
directly to low resolution observations, the obser- 
vations and models were smoothed to a resolving 



power of i?=100 before fitting them. 

• Continuum extinction. The goal of this work is 
to measure the strength and shape of dust and 
ice absorption features, and therefore a continuum- 
only extinction curve should be applied to the 
stellar models. Published extinction curves are 
usually based on broad-band photometry (e.g., 
Ilndebetouw et al.ll2005l: iChapman et al.ll2009[ ). in- 
cluding continuum and feature extinction and ap- 
plying large interpolations between, e.g., the IRAC 
8 /im and MIPS 24 /im photometry. Therefore, 
an independent "high resolution" extinction curve 
is derived f tj4.2l) . from which a continuum-only ex- 
tinction curve was extracted (dashed red line in 
Fig- [2b)- The latter was applied to the stellar mod- 
els with as a free parameter. 
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Table 2 

Model Spectra 



Spectral Type 




log(9) 




K 




G8 III 


5000 


2.50 


KO III 


4500 


2.50 


K3 III 


4500 


2.00 


K4 III 


4250 


2.00 


K5 III 


4000 


2.00 


K7 III 


4000 


1.50 


MO III 


3700 


1.00 


Ml III 


3500 


1.00 


M3 III 


3200 


1.00 


M6 III 


3200 


0.50 


M7 III 


3000 


0.50 


M8 III 


2800 


0.50 


M9 III 


2500 


0.50 


Note. — Spectral typ 


cs were 



assigned using polynomial fits to 
empirical T^ff and \og{g) data fol- 
lowing the method of I Kane" fc Sahul 

• H2O ice column density. Broad II2O ice features 
cover much of the observed background star spectra 
(Fig. [5]). For sources with observed L-band spec- 
tra, the strong, relatively sharp 3.0 /.tm band (O- 
H stretch mode) constrains the H2O ice column 
density well. This significantly improves the relia- 
bility of the overall fits, because in the 5-8 and 11- 
20 /.tm ranges, H2O features overlap with other fea- 
tures. Optical constants of arnorph ous solid II2O 
at T = 10 K fHudgi ns et"all |19 9 3') were used to 
calculate the absorption s pectrum of ice spheres 
(|Bohren fc Huffmaiilll983D . Spheres with radii of 

0. 4 /im fit the typical short wavelength profile and 
peak position of the observed 3 /xm bands best. 
The peak strengths of the 6.0 /im bending and 13 
/im libration modes relative to the 3.0 stretching 
mode of this calculated spectrum are ~7% weaker 
with r espect to the laborat ory transmission spec- 
trum (jHudgins et al.l 119931 ). Also, the libration 
mode peaks at 12.5 /im, compared to 13.3 /im in 
the laboratory. 

• Silicate feature depth. The contribution of sil- 
icate features to the absorption is quantified 
by the peak optical depth of a synthetic sili- 
cate spectrum at 9.7 /im (rg.y). The silicate 
spectrum was calculated for grains small com- 
pared to the wavelen gth (Rayleigh l i mit) us- 
ing optical constants (|Dorschner et al.l 119951 ) of 
amorphous pyroxene (Mgg jFeo.sSiOa) and olivine 
(Mg]^FeiSi04). The olivine and pyroxene spectra 
were added such that the shape of the silicate fea- 
ture in 2MASS J21240517-I-4959100 is matched, 

1. e., Tg.7(pyroxene)/T9.7(olivine)=0.62. Note that 
this silicate spectr um is broader than the diffuse 
medium spectrum (jKemper et al.|[2"004D . for which 
one would need a ratio of 0.20. 

values were derived for all available observables, in- 
cluding photometric points, ground-based spectra and 
Spitzer spectra, except for those for which it is a pri- 
ori known that insufficient information is available: the 



3.1-3.7 /tm spectral region (long-wavelength wing 3.0 /tm 
band), the 5.3-7.1 /tm region (6.0 and 6.85 /im ice bands), 
and the IRAC2 photometric point (4.67 /tm CO ice 
band). The MIPSl photometric point was not fitted ei- 
ther, to avoid too much weight to the longer wavelengths, 
which are not important for the ice bands studied here 
(although, in practice, the final models are in quite good 
agreement with the MIPSl photometry; Fig. [Ij. Be- 
cause there are many more spectroscopic than photo- 
metric data points, the latter would have much larger 
values. Thus the spectroscopic values were divided by 
the number of data points and then the "reduced val- 
ues" of the different observables were averaged to repre- 
sent the goodness of fit over all observables. In addition, 
a reduced value was derived for the most crucial ob- 
servables for spectral type determination: near-infrared 
photometry, the shape of near-infrared spectra and the 
4.0-4.1 /tm spectral region. The reduced values and 
parameters of the best fits are reported in Table [3] and 
the fits are plotted in Fig. [1] (red lines in the left pan- 
els). Good fits are often obtained. In seven cases, the 
reduced x^ values exceed 3 and the main cause for that 
is indicated in the notes of Table [S] Typically, the rela- 
tive scalings of the K-band photometry, and the (partial) 
K-band and L-band spectra deviate from the model, but 
the fits are good at longer wavelengths. 

In seven other cases, no spectra below 5 /tm, and thus 
no independent measures of spectral type and Af(H20) 
were available (Table [iJ. This most strongly affects the 
iV(Il20) determination, as the depth of the 13 /tm li- 
bration mode is degenerate with the slope of the applied 
extinction curve and the assumed silicate band profiles. 
An uncertainty in the spectral type affects the shape and 
depth of the 5-8 nm absorption features because of photo- 
spheric features that emerge for spectral types later than 
Ml III. In practice, the 8 /tm photospheric SiO bands can 
be used to constrain the spectral type sufficiently well, 
however. 

4.2. The "High Resolution" Extinction Curve 

The background star spectra can be used to de- 
rive a high resolution extinction curve, in which the 
ice and dust features can be separated from the con- 
tinuum extinction and in which large interpolations 
usually applied between, e.g., the IRAC 8 and 
MIPS 24 /im photometric data points are avoided. 
2MASS J21240517-h4959100 behind the core L1014 was 
chosen for this purpose, because high quality near- 
infrared spectra are available, indicating an accurate 
spectral type of Ml III. The extinction curve was de- 
rived as follows: 

AA[mag] = -2.5Log(i^A(obs)/FA (model)). (1) 

The absolute scaling of F\ (model) is unknown, because 
there are no measurements of the star's distance. To 
circumvent this problem, the model was scaled with re- 
spect to the observed spectrum such that the derived 
Ax follows known near- infrared extinction curves (i.e.. 
Ax oc A-i ®; e.g.. Ilndebetouw et al.ll2005D . This pro- 
vides the "high-resolution" extinction curve shown in 
Fig. [2]d. This curve is remarkably fiat. At an Ak of 3.02 
mag, the extinction is still 1.5 mag at 25 /tm. This con- 
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Figure 1. Left panels: Observed ground- based and Spitzer/IRS spectra combined with broad band photom etry (filled circles), and lower 
limits (open triangles) and 3(t upper limits (open circles) thereof. The red line represents the fitted model ( i|4.ip . The sources are sorted 
in increasing Right Ascension. The name of the background stars and the cores they trace are indicated. Right panels: Optical depth 
spectra, derived using the continuum models of the left panels, excluding the modeled ice and silicate features. Note the large variation of 
the relative depths of the 3.0 and 9.7 /im bands. The red line indicates the modeled H2O ice and silicates spectrum. For clarity, error bars 
of the spectral data points are not shown. 
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Table 3 

Continuum Fit Parameters 



Source 


Spectral Type" 






T9.7'' 


X?(total)° xUSpTY 


2MASS J 




mag 









r\A 01 /I no 


1 1 t^QHonn 


KO (G8-K5) 


Q 


HA 
U4 




uy ) 


1 no fn nf^A 
l.Uy (^U.UoJ 


U.oo 


/n n/i ^ 
(^U.U4 ) 


1 on 


i 


y3 


f\Ci r\K or* 

08052135 


— 3909304 


MO 


(K5-Ml)e 


1 


47 


(0 


04 ) 


<0.70'' 


0.41 


(0.02 ) 


0.29 






08093135 


- 3604035 


KO 


'G8-Ml)s 


1 


89 


(0 


06 ) 


0.55'' (0.25) 


0.58 


(0.02 ) 


0.42 






08093468 


- 3605266 


M3 


M0-M7)s 


4 


52 


(0 


14 ) 


2.10'' (0.42) 


0.95 


(0.04 ) 


0.58 






12014301 


- 6508422 


M3 


M3-M7)s 


5 


17 


(0 


15 ) 


3.201^ (0.80) 


1.35 


(0.05 ) 


1.94 






12014598 


- 6508586 


K3 


(G8-K7)s 


3 


07 


(0 


09 ) 


1.45'^ (0.25) 


0.78 


(0.03 ) 


0.37 






15421547 


- 5248146 


MO 


(K7-Ml)g 


1 


94 


(0 


06 ) 


<0.70'^ 


0.81 


(0.03 ) 


0.51 






15421699 


- 5247439 


K7 


K5-M3)s 


2 


87 


(0 


09 ) 


0.85'^ (0.25) 


0.84' 


(0.03 ) 


0.26 






17111501 


- 2726180 


M6 


(M3-M6) 


3 


91 


(0 


12 ) 


1.39 (0.07) 


0.76 


(0.03 ) 


1.21 


1 


92 


17111538 


- 2727144 


M3 


{M1-M6) 


3 


53 


(0 


11 ) 


1.20 (0.06) 


0.87 


(0.03 ) 


1.96 


1 


69 


17112005 


- 2727131 


M3 


(M1-M3) 


6 


00 


(0 


18 ) 


2.25 (0.11) 


1.30 


(0.05 ) 


4.91'' 


5 


22 


17155573 


- 2055312 


Ml 


(MO-Ml) 


2 


65 


(0 


08 ) 


0.97 (0.05) 


0.66 


(0.03 ) 


1.57 


1 


94 


17160467 


- 2057072 


K7 


(K7-M0) 


3 


24 


(0 


10 ) 


1.27 (0.06) 


0.77 


(0.03 ) 


9.38'' 


3 


42 


17160860 


- 2058142 


G8 


(G8-K4) 


2 


45 


(0 


07 ) 


0.91 (0.05) 


0.77 


(0.08 ) 


0.72 


1 


28 


18140712 


- 0708413 


K7 


{K5-M0) 


1 


89 


(0 


06 ) 


0.83 (0.04) 


0.60 


(0.02 ) 


0.64 





72 


18160600 


- 0225539 


MO 


(MO-Ml) 


1 


25 


(0 


04 ) 


0.63 (0.09) 


0.50 


(0.02 ) 


1.77 


1 


76 


18165296 


— 1801287 




Ml 


3 


00 


(0 


18 ) 


0.96 (0.10) 


2.00 


(0.12 ) 


1.21 


1 


34 


18165917 


- 1801158 


M6 


{M3-M6) 


3 


34 


(0 


10 ) 


1.10 (0.05) 


1.66 


(0.07 ) 


0.76 


1 


10 


18170426 


- 1802408 


M6 


(M3-M6) 


2 


50 


(0 


08 ) 


0.65 (0.03) 


1.38 


(0.06 ) 


1.89 


4 


22 


18170429 


- 1802540 


Ml 


(MO-Ml) 


2 


96 


(0 


09 ) 


0.85 (0.09) 


1.40 


(0.06 ) 


2.07 


1 


92 


18170470 


- 0814495 


MO 


(MO-Ml) 


4 


27 


(0 


13 ) 


2.33 (0.12) 


1.39 


(0.06 ) 


12.28'^ 


2 


49 


18170957 


- 0814136 


Ml 


(M1-M6) 


3 


27 


(0 


10 ) 


1.69 (0.08) 


1.04 


(0.04 ) 


2.02 


2 


49 


18171181 


~ 0814012 


K7 


(K5-M0) 


4 


31 


(0 


13 ) 


2.26 (0.11) 


1.49 


(0.06 ) 


2.00 


1 


11 


18171366 


~ 0813188 


K7 


(K5-M0) 


3 


58 


(0 


11 ) 


2.02 (0.10) 


1.17 


(0.05 ) 


16.21"^ 


1 


93 


18172690 


~ 0438406 


M3 


(M1-M6) 


4 


60 


(0 


14 ) 


2.55 (0.13) 


1.11 


(0.04 ) 


4.53'^ 





78 


19201597 


+ 1135146 


Ml 


(MO-Ml) 


3 


14 


(0 


09 ) 


0.79 (0.04) 


1.05 


(0.04 ) 


1.01 





69 


19201622 


+ 1136292 


M6 


(M3-M6) 


2 


50 


(0 


13 ) 


0.60 (0.03) 


0.78 


(0.03 ) 


0.99 


1 


11 


19214480 


+ 1121203 


M6 


(M3-M6) 


2 


10 


(0 


06 ) 


0.85 (0.04) 


0.95 


(0.04 ) 


7.31'^ 


1 


74 


21240517 


+ 4959100 


Ml 


(MO-Ml) 


3 


10 


(0 


09 ) 


1.30 (0.06) 


0.95 


(0.04 ) 


1.03 





82 


21240614 


+ 4958310 


K4 


(K4-K5) 


1 


60 


(0 


05 ) 


0.58 (0.03) 


0.52 


(0.02 ) 


1.35 


1 


33 


22063773 


+ 5904520 




G8 


1 


74 


(0 


09 ) 


0.67 (0.03) 


0.33 


(0.01 ) 


2.95 


2 


48 


04393886 


+ 2611266 


K7 


(K5-M1) 


3 


00 


(0 


09 ) 


1.42 (0.07) 


0.84 


(0.03 ) 


0.80 


2 


61 


18300061 


+ 0115201 


M3 


(M3-M6) 


6 


38 


(0 


19 ) 


1.81 (0.09) 


1.71 


(0.07 ) 


7.30'^ 


18.79 



^ Best fitting spectral type that is used throughout this work. The uncertainty range is given in parentheses, unless 
the spectral type is more accurate than the closest available sub-types. Note that the spectral types are limited to the 
ones listed in Table [2] All listed spectral types have luminosity class III. 
^ Extinction in the /C-band 

^ Peak absorption optical depth of the 3.0 f-Lm. H2O ice band. 
^ Peak absorption optical depth of the 9.7 band of silicates. 

^ Reduced values of the model spectrum fitted to all available photometry and spectra. 

^ Reduced values of the model spectrum to all available near-infrared photometry and spectra (J, H, K, and 
L-band), excluding the 3.0 fim ice band. 

^ Spectral type determined from wavelengths above 5 fim only because no near-infrared spectra available. 

^ Ta.o uncertain because no near-infrared spectra available. 

^ 9.7 fim band narrower than for most other sources in the sample. 

J No near-infrared spectra available. 

^ The main reason for large values of x^(total) is different in individual cases: 2MASS J17112005-2727131: only partial 
K-band spectrum available that does not match K-band photometry, and only upper limits to 2MASS J and H-band 
photometry available. 2MASS J17160467-2057072: only partial K-band spectrum available that does not match K-band 
photometry, possibly due to flux calibration error of the spectrum. 2MASS J181704 70-0814495: combined partial K and 
L-band spectra too high compared to model, uncertainty in IRACl photometry underestimated? 2MASS J18171366- 
0813188: K-band spectrum does not match K-band photometry, possibly relative K versus L-band calibration error. 
2MASS J1S172690-043S406: model too shallow between the H and K-bands, and too steep above 10 fim. 2MASS 
J192144S0-H1121203: modeled K- and L-band spectra systematically too low. 2MASS JlS300061-t-0115201: no IRAC 
photometry available for scaling the spectra, H and K-band slope steeper than can be modeled. 



firms the results of lChapman et al] (I2Q0S ), as in dicated 
by the symbols in Fig. Eb. and IMcClur j ((1009). Sub- 
sequently, continuum and feature extinction were sepa- 
rated by ad ding; a laboratory spe ctrum of soHd H2O (at 
T = 10 K; iHudgins et al.l I1993D and a siHcates model 
spectrum ( N4.ip to a polynomial fit to feature-free re- 
gions of the form 



Log{Ax/AK) - ao + aiLog(A) -h a2[Log(A)]^ 



(2) 



with A in /im and Ax in magnitude. The polynomial 
coefficients are listed in Table HI and the curve is plotted 
with a dashed red line in Fig. [2b. This feature-free ex- 
tinction curve works well for most sources. Examples of 
exceptions might be 2MASS J19201597-I-1135146 (model 
too shallow > 10 Atm) and 2MASS J18172690-0438406 
(model too steep > 10 /im), but no systematic trend with 
or core environment could be found. In these cases 
the silicate model used in the fits may be responsible for 
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Figure 2. Derivation of an extinction curve for the star 
2MASS J21240517-I-4959100 behind L1014. Panel (a) shows the 
observed spectrum, broad-band photometry (diamonds) and an 
Ml III model spectrum. After subtracting these spectra on a Log 
scale, the resulting extinction curve is shown in panel (b). The 
solid red line is the result of adding spectra of solid H2O (thick solid 
black line) and silicate features (dotted black line) to a polynomial 
fit. The dashed line then represents the feature-free extinction 
curve used in the background star continuum fits. Green sym- 
bols represent the average extin ction in broad band fil ters derived 
from nearby cloud observations l|ChaDman et al1l2009i ) for lines of 
sight with < 0.5 (open squares), 0.5 < j4k < 1-0 (open cir- 
cles), 1.0 < tIk < 2.0 (filled squares), and 1.0 < < 2.0 (filled 
circles). Panel (c) compares the empirical extinction curve with 
the calculated curves of /?v= 3.1 (case A) and 5.5 (case B) from 
IWeingartner Drain^ 1I2OOII '). 



Table 4 

Polynomial 
Coefficients of 
Feature- Free 
Extinction Curve 
(Eq.[3 



Coefficient Value 



ao 


0.5924 


ai 


-1.8235 


a2 


-1.3020 


as 


5.9936 


a4 


-5.3429 


as 


1.2619 


ae 


0.2738 


ay 


0.0069 


ag 


-0.0554 



the deviations as well. 

4.3. Absorption Band Strengths and Column Densities 

The continuum fits, as indicated with red lines in Fig.[T] 
(left panels), but excluding modeled H2O and silicate 



bands, were used to put the spectra on an optical depth 
scale. The right panels of Fig. [1] show that the optical 
depth in the 7.5-8.0 /im region is not always fully ex- 
plained by the modeled H2O and silicate bands. This 
may be related to inaccuracies in the applied extinction 
curve and stellar model. Therefore, in the subsequent 
analysis of the optical depth spectra, local straight line 
continua between 7.5 and 5.3 or 6.4 /im were used to 
correct for these small offsets (see Fig. [3] for examples). 

Peak and integrated optical depths were determined 
for the 6.0 and 6.85 /im bands, before and after the 0-H 
bending mode of H2O was sub tracted using an am or- 
phous ice spectrum at T^IO K (|Hudgins et al.|[T993l : Ta- 
ble [5]). The features in the H20-subtracted spectra were 
also decomposed int o the five components discussed in 
iBoogert etall (|2008[ ). These components (labeled "Cl"- 
"C5" ) are useful in characterizing source-to-source vari- 
ations of the absorption in the 5-7 /im wavelength region 
( ij4.4.6p . They do not represent a unique decomposition, 
however, and each component may be the result of ab- 
sorption by more than one molecular band. Examples 
of this decomposition are shown in Fig. |31 For compar- 
ison two YSOs are shown as well (Boogert et al. 2008). 
Clearly, the relative depths of the components vary con- 
siderably between background stars and between YSOs. 

Column densities were determined for solid H2O, NH^, 
and CO2 (Table lU, using the intrinsic integrated band 
strengths listed in Table [7l The II2O column density 
is unreliable if no L— band spectrum is available (7 out 
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Figure 3. Decomposition of the absorption features in the 5-8 
/im spectra of two background stars (bottom panels) , compared to 
two YSOs ( Boogert et al. 2008) in the top panels. The components 
are labeled CI (solid green), C2 (dashed green), C3 (solid purple), 
C4 (dashed purple), and C5 (dashed red). The sum of the com- 
ponents is indicated with a thick red line. Note that the relative 
strengths of the components (e.g., C3/C1) varies strongly toward 
both YSOs and background stars. Component C5 has not been 
detected toward any background star. All components are added 
to a local straight line continuum (dotted black line) to remove 
small inaccuracies from the global continuum fitting process. 



Table 5 

Optical Depths 5-8 fim Features [publish electronically] 
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Tint [cm 1] 
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Note. — Uncertainties in parentheses based on statistical errors in the spectra only, unless noted otherwise below, 
integrated optical depth between 5.2-6.4 /im in wavenumber units 
^ integrated optical depth between 5.2-6.4 /jm in wavenumber units, after subtraction of a laboratory spectrum of pure H2O ice 
integrated optical depth between 6.4-7.2 fira in wavenumber units 

integrated optical depth between 6.4-7.2 fim in wavenumber units, after subtraction of a laboratory spectrum of pure H2O ice 

peak optical depth at 6.0 fira 

' peak optical depth at 6.85 /im 

s peak optical depth component CI IIBoogert et al.ir200SI1 
^ peak optical depth component C2 
' peak optical depth component C3 
J peak optical depth component C4 

^ peak optical depth component C5. Uncertainty includes uncertainty in A'^(H20) as this affects the baseline level. 
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2 


41 


(0.42) 


6.13 (1.26) 


2.49 


(0.60) 


6.34 (1.69) 


<1.28 


<3.66 


<1 


80 


<5.14 


<7.00 


<20.00 


18170957 


- 0814136 


2.85 


(0.31) 


1 


82 


(0.27) 


6.39 (1.20) 


12.29 (1.29) 


43.12 (6.63) 


2.69 


(0.58) 


9.46 (2.30) 


<0.89 


<3.53 


<0 


97 


<3.84 


<3.16 


<12.48 


18171181 


- 0814012 


3.81 


(0.42) 


2 


80 


(0.50) 


7.35 (1.55) 






4.42 


(0.86) 


11.58 (2.60) 


<1.48 


<4.38 


<1 


42 


<4.19 


<6.71 


<19.80 


18171366 


- 0813188 


3.40 


(0.38) 


3 


09 


(0.48) 


9.08 (1.74) 






3.48 


(0.86) 


10.24 (2.78) 


<0.99 


<3.29 


<1 


24 


<4.12 


<4.83 


<15.97 


18172690 


- 0438406 


4.31 


(0.48) 


3 


02 


(0.30) 


7.01 (1.05) 


18.86 (1.72) 


43.75 (6.31) 


3.63 


(0.65) 


8.44 (1.79) 


<1.63 


<4.28 


<0 


95 


<2.48 


<7.00 


<18.28 


19201597 


+ 1135146 


1.34 


(0.14) 


1 


42 


(0.39) 


10.64 (3.14) 








<0.70 


<5.95 


<0.37 


<3.12 


<1 


15 


<9.69 


<6.58 


<55.33 


19201622 


+ 1136292 


1.01 


(0.11) 





31 


(0.17) 


3.09 (1.74) 


3.86 (1.07) 


38.08 (11.3) 




<0.70 


<7.77 


<0.50 


<5.58 


<0 


58 


<6.49 


<2.00 


<22.15 


19214480 


+ 1121203 


1.43 


(0.16) 


1 


38 


(0.50) 


9.65 (3.68) 




<1.61 


<12.64 


<0.49 


<3.88 


<1 


63 


<12.85 


<2.57 


<20.22 


21240517 


+ 4959100 


2.19 


(0.24) 


1 


48 


(0.29) 


6.76 (1.55) 


8.24 (2.69) 


37.60 (13.0) 


1.54 


(0.38) 


7.04 (1.90) 


<0.40 


<2.07 


<1 


83 


<9.43 


<4.20 


<21.57 


21240614 


+ 4958310 


0.98 


(0.11) 


1 


18 


(0.34) 


12.01 (3.78) 








<0.59 


<6.83 


<0.27 


<3.13 


<1 


14 


<13.04 


<2.50 


<28.57 


22063773 


+ 5904520 


1.13 


(0.12) 


1 


44 


(0.39) 


12.70 (3.74) 








<0.40 


<3.98 


<0.31 


<3.10 


<1 


11 


<11.02 


<2.30 


<22.82 


04393886 


+ 2611266 


2.39 


(0.26) 


1 


44 


(0.06) 


6.03 (0.72) 


7.11 (0.05) 


29.66 (3.32) 




<0.25 


<1.17 


<0.38 


<1.81 


<0 


82 


<3.84 


<1.45 


<6.85 


18300061 


+ 0115201 


3.04 


(0.34) 


2 


68 


(0.20) 


8.79 (1.18) 


11.93 (0.21) 


39.14 (4.43) 




<2.10 


<7.75 


<1.18 


<4.36 


<0 


58 


<2.14 


<6.39 


<23.61 



a 

fD 




o 

o 



Note. — Column densities were determined using the intrinsic integrated band strengths listed in Table [7] Uncertainties (Itr) are indicated in brackets and upper limits are 
of 3(T significance. 

^ An uncertainty of 10% in the intrinsic integrated band strength is taken into account in the listed column density uncertainties. 

^ Assuming both C3 and C4 components lIBoogert et al.ir2008l ) are due to NHJ ( i|4.3|l . The contribution by the CH3OH C-H bending mode has been subtracted for sources 
with CH3OH detections. 

The uncertainties in the CH3OH column densities do not reflect interference with photospheric absorption lines (Fig.|4]|. CH3OH detections toward 2MASS J18171181-0814012 
and 2MASS J18171366-0813188 are most reliable. 

No L-band spectrum available for H2O ice column density determination. The 13 /im H2O libration mode was used instead. These values arc uncertain, and the listed 
uncertainties are estimates. 
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Table 7 

Intrinsic Integrated Band Strengths 



Molecule 


Mode 


A 


Reference 






cm / molecule 





H2O 3.0 /im stretch 

NH4' 6.8 /xm deformation 

CO2 15.0 lira bend 

CH3OH 3.53 /tm C-H stretch 

CH3OH 9.7 /im O-H stretch 

HCOOH 5.85 lira C=0 stretch 

HCOOH 7.25 lira C-H deformation 

CH4 7.68 lira deformation 

NH3 8.9 lira umbrella 



2.0 X IQ-lS 
4.4 X 10"" 

1.1 X 10"" 
5.6 X 10-18 

1.6 X 10"" 

6.7 X 10"" 
2.6 X 10-1* 
7.3 X 10-1* 
1.3 X 10-" 



Hagen et al. (1981 1 
Schutte fc K hanna (2003) 
Gerakines et al. (1995) 



Kerkhof et al. (1999) 
Kcrkhof et al. (1999) 
Schutte et ah (1999 ) 
Schutte et ar (1999 ) 
Boogert et aL (1997 ) 
Kerkhof et al. (1999) 




10 
(/xm) 

Figure 4. Left panel: Sources with detections (top 2 sources) or tentative detections of the C-H stretching mode of solid CH3OH 
at 3.53 lira . The green line is t he assumed local continuum and the red line the laboratory spectrum H2O:CH3OH:CO:NH3=100:10:l:l 
(T = 10 K: IHudgins et anil993l ) corresponding to the column densities listed in Table [6l Middle panel: the 6.3-8.0 lira region of the same 
sources as in the left panel showing the C-H bending mode of CH3OH (red) with the assumed local continuum (green) at the same column 
densities as in the left panel. The observed absorption is much deeper due to an overlapping band of (possibly) NH^ ( ^4.3| l. Right panel: 
the 8-11.5 ^m region of the same sources as in the left panel showing the C-O stretching mode of CH3OH (red) with the assumed local 
continuum (green) at the same column densities as in the left panel. Note that a scaling factor along the optical depth axis was applied 
for the spectra to fit in the panel. 



of 31 sources), and the 13 /.tm libration mode is the 
main tracer. The depth of the latter is strongly af- 
fected by the assumed shape of the extinction curve 
and by overlapping stretching and bending modes of 
silicates. The NH4 column density is derived, assum- 
ing that the entire 6.85 fim. band (i.e., both compo- 
nents "C3" and "C4" listed in Table O can be at- 
tributed to it, after correction for any overlapping bands 
of the C-H deformation mode of solid CH3OH (see be- 
low). This identification has not been fully validated, 
howeve r, because no other NHj" bands have been de - 
tected (jSchutte fc Khannal 120031 : iBoogert et all I2008D . 



The identification of the 15.2 /im band is not disputed, 
and the solid CO2 column density is listed for the 8 
sources for which this wavelength range was observed. 
All CO2 observations were done at a low resolving power 
of ~90, which is insufficient to assess the composition or 
thermal processing history of the ices (jPontoppidan et aTl 
[2OO8I) . 

The background star spectra were also searched for 
CH3OH, HCOOH, NH3 and CH4 species that are often 
detected toward YSOs. For solid CH3OH, the 3.53 ^m 
C-H stretching and the 9.7 /xm C-O stretching modes 
are usually used for column density determinations 
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(|Boogert et al.|[2008l:[Bottinelli et alJlMol) . The 3.53 /im 
band is preferred because it is free from interfering strong 
bands. Neverthel ess it overlaps with the 3.47 /im band 
of hy drocarbons (jAllam andola et al. 1992; Brook e et al.l 
I1996D or NH3 hydrates (partois fc d'Hendecourtll200lD . 
as well as several photospheric lines. Thus, the most cer- 
tain detections of solid CII3OII are toward background 
stars that show a distinct 3.53 /im band. Fig. |4] shows 
that this is the case for certainly two lines of sight: 
2MASS J18171181-0814012 and 2MASS J18171366- 
0813188. The strongest detection (2MASS J18171181- 
0814012) was verified by observations in two different ob- 
serving runs, using two different telluric standard stars. 
A good agreement is found for the strength of the 9.7 /xm 
band. Five additional lines of sight likely show solid 
CII3OH as well, although in these cases the 3.53 /im 
feature is not as distinct from other fluctuations in the 
surrounding spectra (Fig. \^ . For all of these 7 sources 
the CII3OII abundance relative to H2O is in the 5-12% 
range. A similar number of sources have upper limits of 
comparable magnitude, and thus the CII3OII abundance 
varies in different environments. 

For the YSOs studied in IBoogert et all (|2008[ ). 
HCOOH is considered a detection only if the 7.25 /im 
C-H deformation mode has been detected. This is not 
the case for any of the background stars. Based on 
the strength of the CI component (Table [5]), which may 
mostly be due to the C=0 stretch mode of HCOOH, the 
abundance can in many lines of sight be constrained to 
less than or comparable to the typical detections toward 
YSOs of 2-5% relative to H2O, however (Tabled). 

Finally, for NH3 and CH4 the 8.9 /j,m umbrella and 
7.68 /im bending modes are used to determine column 
density upper limits. For CH4 the upper limits are some- 
times comparable but rarely less than the detections of 
-4% toward YSOs (|Oberg et al.ll200l . For NH3 the up- 
per limits are typically ~20% relative to H2O (Table 
and are thus not signifi cant compared to the detections 
of 2-15% toward YSOs (jBottinelli et al.llMol ). 

4.4. Correlation Plots 

To facilitate the analysis of the observational param- 
eters in many sight-lines in many different environ- 
ments (quiescent core and cloud media, low mass YSOs, 
massive YSOs), a number of correlation plots are pre- 
sented below. The v alues for the YSOs were taken from 
IBoogert et aH (|2008[ ) and lReach et al.l (|2009( ). Note that 
in plots involving Ak no YSOs are present because of 
the unknown shape of the intrinsic K-band continuum 
emission. Also, background stars with T3.0 upper limits 
have been excluded from the plots. 

4.4.1. Tg.r versus Ak 

The plot of the peak optical depth of the 9.7 fj,m 
band of silicates against K-band continuum extinction 
Ak shows an interesting behavior ( Fig. [5^). Wh ile at 
low Ak the diffuse medium relation ()Whittetl 120031) 

r9.7 = 0.554 X Ak (3) 

is followed (assuming Ay /Ak=8.8), a less steep relation 
is observed at Ak >1 mag in dense sight-lines. This was 



discovered bv lChiar et all (|200l in the IC 5146, B68 and 
B59 cores and Serpens and Taurus clouds, and here this 
behavior is extended in more lines of sight, at higher K- 
band extinctions. Some background stars, however, fol- 
low the diffuse medium correlat ion, most no tably those 
behind L328, and, as noted bv Chiar et al.l ([2007), one 
source behind Serpens (2MASS J18285266-H0028242). 
Thus, in these lines of sight the extinction may well be 
dominated by diffuse rather than dense dust. Exclud- 
ing these particular background stars, a polynomial was 
fitted to objects with Ak >1.4 mag: 

rg.T = (0.12 ± 0.05) + (0.21 ± 0.02) x Ak (4) 

It is worth noting that for some cores all background 
stars lie systematically above or below the fit, e.g., L429- 
c sources lie above it and B59 sources below it. 

4.4.2. ra.o versus Ak 

The peak optical depth of the 3.0 /im H2O ice band 
correlates strongly with Ak (Fig. [SJj). There is no indi- 
cation of flattening. A linear fit to all data points yields: 

T3.0 = (-0.15 ± 0.13) 4- (0.45 ± 0.05) x Ak (5) 

This relation implies a T3.0 — cut-off value of Ak — 
0.33 ± 0.30, which corresponds to A\r = 2.9 ± 2.6 (as- 
suming Ay /Ak=8.8). This is the so-called ice formation 
threshold, and it is comparable to the value observed for 
the T aurus Molecular Cloud {Ay = 3.1±0.6: [CIiiar et al.l 
Il995f). but smaller tha n the Ay = 10 — 15 mag quoted 
bv lTanaka et al.l (|1990l ) for Ophiuchus. The large uncer- 
tainty must refiect different ice formation thresholds in 
the different environments traced by the observations, or 
the presence of different amounts of ice-less diffuse dust 
along the line of sight. For example, all L328, B59 and 
CB 188 sources lie to the right of the fitted line. For 
L328 this most likely indicates a significant contribution 
of extinction in diffuse foreground clouds as indicated by 
the Tg.y versus Ak relation f ii4.4.ip . On the other hand, 
all L429-C sources lie to the left of the fit, indicating a 
low ice formation threshold. 

4.4.3. ra.o versus rg.r 

Although the general relation of ts.q with T9.7 is evident 
from panels (a) and (b) of Fig. \5\ it is plotted separately 
in panel (c). A linear fit to all points except those from 
L328 f !HXT|) yields 

T9.7 = (0.36 ± 0.09) + (0.36 ± 0.06) x T3.0 (6) 

For this p a.rticular relation, valu es for the YSOs are avail- 
able from IBoogert et all (|2008D . They are displayed in 
Fig. ini where the linear fit from Eq. [S] above has been 
plotted. YSOs deviating more than 3cr from the fit are 
labeled. YSOs lying below the line may have some of the 
silicate band filled in by emission, while sources above 
the line may have significant amounts of warm dust in 
which the H2O ices have evaporated. Interestingly, the 
latter are mostly the massive YSOs and YSOs that show 
a "red" 6.85 /im band, indicative of thermally processed 
ices (gXg). 
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Figure 5. Correlation plots of with the peak optical depth of the 9.7 /im feature (a), of with the peak optical depth of the 3.0 fim 
feature (b), and of the peak optical depths of the 3.0 and 9.7 fim features (c). Background stars tracing different clouds and cores are 
indicated with different symbols and colors. Solid black lines show least-square fits to all data points, excluding the ones strongly affected 
by diffuse medium absorption (L328) and the ones in the Taurus and Serpens molecular clouds. The fit in panel (a ) was only done for 
Ak > 1.4 mag. The dashed line in that panel is the diffuse medium relation taken fro m the literature llWhitteti 12003) . and the values for 
the B68, IC 5146, Chamaeleon, and most Serpens and Taurus sources were taken from lChiar et aEI 120071) . 
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Figure 6. Peak optical depths of the 3.0 and 9.7 fj,m absorption features for both background stars (in black, red and green; see legend 
of Fig. [Sjl and YSOs (purple asterisks). Only YSOs with reliable 3 fim spectra are plotted. The solid line represents the fit for background 
stars as derived in Fig. [5]panel (c). YSOs deviating more than 3cr from the solid line are labeled, except for IC 1396 a and B5 IRSl (both 
above the solid line), and IRAS 03439, L1448 IRSl (below the solid line), which are not labeled to avoid crowding in the plot. 
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Figure 7. Integrated optical deptlis of H2 0-subtractcd spectra 
in the 5.2-6.4 (panel a) and 6.4-7.2 (panel b) ranges plotted 
against the H2O ice column densities. The symbol types and colors 
have the same meaning as in Figs. \E\ and |6] Despite the large 
scatter, it can be concluded that the background stars follow the 
same correlations as the YSOs (purple asterisks), i.e., the carriers 
of the absorption are already present in quiescent cloud material. 

4.4.4. 6.0 fim band versus A'^(H20) 



Prev ious work by Keane et al.l ()2001bD toward massive 
YSOs. iK^z et all (I2005D toward background stars and 
IBoogert et alj (|2008[ ) toward low mass YSOs has shown 
that the 6.0 /im band can only partly be attributed to 
solid II2O. Although H2O mixed with realistic concen- 
trations of C O2 fits the long wa velength side (compo- 
nent "C2" in IBoogert et"alll2008l) be tter than pure H2O 
(jKnez et al.ll2005HOberg et al.ll2007| ). it does not explain 
all absorption. This is the case for the background stars 
studied here as well. To quantify this, the integrated ab- 
sorption remaining after subtraction of an amorphous, 
pure H2O laboratory ice at T =10 K (jHudgins et al.l 
I1993D is plotted against A^(H20), showing that the back- 
ground stars are in line with the YSOs (Fig.[7K; Tabled]). 
The scatter is large, however, which is refiected in a lin- 
ear correlation coefficient of 0.79. Some of this might 
be due to uncertainties in A^(H20), as many YSOs do 
not have reliable 3.0 /im spectra available. Indeed, when 
only including sources with reliable 3.0 /um spectra, the 
correlation coefficient increases to 0.85. Another source 
of scatter may be the inaccuracy in the continuum deter- 
mination, especially for the YSOs. Intrinsic abundance 
variations of the different carriers of the 5.2-6.4 fim "ex- 
cess" absorption relative to II2O must cause some of the 
scatter as well (Fig. [S]). 

4.4.5. 6.85 fim band versus A'^(H20) 

The 6.85 /im band is detected toward many back- 
ground stars. Its integrated optical depth, calculated 
after subtracting the long wavelength wing of the II2O 
bending mode, is correlated with A^(H20) in Fig. [7|d, 
where the YSOs have been plotted as well. As for the 
6.0 /im band, the scatter is large. At a given value of 
A^(H20), the 6.85 /xm band depth varies by a factor of 
2. The linear correlation coefficient is 0.88 when includ- 
ing all background stars and YSOs. Within the scatter, 
there is no systematic difference in the strength of the 




T. ,C4/t. ,C3 
mt ' mt 

Figure 8. Ratio of integrated optical depths of the C4 and C3 
components (long and short wavelength components of the 6.85 fim 
band) plotted versus the H2O ice column density normalized to the 
peak optical depth of the 9.7 ^m band (a measure of the H2O ice 
abundance). Only values with significance > 2cr are plotted. Dif- 
ferent colors and symbols refer to background stars behind different 
clouds and cores (see legend of Fig.fSjl. The valu es for the YS Os, in- 
dicated with purple asterisks, are taken from iBoogert et al.l f2008i ) . 
While nearly all sources have C4/C3 ratios near 1, a number of 
YSOs have significantly larger ratios and at the same time low 
H2O abundances. 

6.85 /im band between YSOs and background stars. 
4.4.6. C4/C3 Ratio versus H2O abundance 



The 6.0 and 6.85 um bands were decom posed into com- 
ponents C1-C5 in IBoogert et al.l (|2008f ). It was found 
that the C4/C3 component ratio varies significantly be- 
tween YSOs. Sources with large C4/C3 ratios tend 
to have low H2O ice abundances (represented by the 
A^(H20)/r9 7 ratio) a nd large deg,r ees of ice crystalliza- 
' all 1200 Ibt IBoogert et al.l 



tion and segregation (jKeane et 

120081) ■ thus indicating that the C4/C3 ratio increases as 
a result of thermal processing. A similar decomposition 
was done for the background stars (Fig. jS]), and indeed 
they have low C4/C3 ratios as expected for the cold, un- 
processed sight-lines that they trace. Nevertheless, some 
of the background stars (e.g., toward L328) have low H2O 
ice abundances, likely because of the presence of diffuse, 
ice-free, foreground clouds ft i4.4.ip rather than because 
of the effects of ice sublimation as for YSOs. 

5. DISCUSSION 



5.1. Ice Abundances 

The infrared spectra of ices and dust toward back- 
ground stars presented here provide a baseline for the 
astrochemical evolution of YSOs. It is shown that the 
strength and shape of the 3.0, 6.0, 6.85 and 15 /im 
ice bands are similar for the quiescent material in iso- 
lated cores traced by background stars compared to most 
previously studied YSOs. Correspondingly, the rela- 
tive abundances of the ices with secure (CO2, CII3OII) 
and less secure identifications (NH^) are similar to most 
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YSOs. Species responsible for the 6.0 band, be- 
sides H2O (probably HCOOH, H2CO, and NH3), must 
have abundances similar to YSOs as well. This con- 
firms earlier work on smaller samples (jKnez et al.l 120051 : 
IBergin et al.l I200"5l) . and it firmly establishes that most 
ices are formed early on in the molecular cloud evolu- 
tion, and are largely unaltered during the further evolu- 
tion of the region (e.g., the star formation process) until 
the YSO thermally processes the ices in the envelope. 

The two background stars with the most secure 
CH3OH detections (2MASS J18171181-0814012 and 
2MASS J18171366-0813188), at ~11% with respect to 
H2O, are both located behind the core L429-C. These 
lines of sight show some of the deepest ice and dust 
features in the sample, which may have facilitated the 
detection of the weak 3.53 /xm CH3OH feature. In- 
deed, many other background stars in the sample have 
upper limits or tentative detections at abundances not 
much less (6-10%) than the detections. Some, e.g. 
2MASS J04393886-h2611266 in the Tauru s Molecular 
Clou d, have tighter upper limits of 1-3% (jChiar et al.l 
119961 : Table[6]), however. Large CH3OH abundance varia- 
tions are also observed toward YSOs (Dar tois ct al. 1999; 
iPontoppidan et al.l[2003n iBoogert et al.lf2008l ). 

The presence of solid CH3OH in quiescent clouds 
is not unexpected, f ollowi n g theoretical studies 
(jTielens fc AUa mandola' 198?; 'Keane '2001a'), labo- 
rator y studies (Watanabc & Kouchi 2002; Fuch s et al.l 
l2009fl and time-dep endent (Monte Carlo) simulations 
(jCuppen et al.l 12009.) . While the laboratory experi- 
ments have not yet been fully translated into kinetic 
parameters, in a general sense, the theoretical studies 
and kinetic Monte Carlo simulations do show that 
the CH3OH/CO ratio is very sensitive to the dust 
temperature a nd to the ratio of atomic H to CO in the 
accreting gas (jKeand l2001al : iCuppen etHI l2009f ). In 
this model, an H20-rich mantle forms first, because of 
its higher sublimation temperature, and at T <17 K a 
CO mantle forms on top of that. Accreting H atoms 
then form H2CO and CH3OH. High CH3OH/CO ratios 
are indicative of high H/CO ratios of the accreting gas 
(due to relatively low densities or high local cosmic ray 
fluxes that destroy H2). Alternatively, low CH3OII/CO 
ratios may r eflect relatively eleva ted dust temperatures 
(T > 15 K: ICuppen eFall [20091 ). Thus local physical 
conditions are a dominating factor in the formation 
of CH3OII, and this may explain the large observed 
CH3OH abundance variations in different lines of sight. 

5.2. Processing of the Ices 

Ice evolution is observed in the envelopes of some 
YSOs, and it appears that this is mostly due to the ef- 
fects of heating: apol ar CO-rich ices sub l imate, and H2O- 
rich ic es crystallize (Smith et al.iri989l: iGerakines et al.l 



1999'; 'Boogert et all 120001: IPontoppidan et al.l 120081 

Bpogcrt ct al. 2008). The carrier of the 6.85 jjim band is 
affected by heating as w ell: it shif ts to longer wavelengths 
(jKeane et a l. 2001b; B oogert et a l. 2008). Toward back- 
ground stars, none of these effects are observed. Fig. [8] 
shows that all background stars with high quality spectra 
have 6.85 /im bands peaking at short wavelengths (low 
C4/C3 ratios). The observed 3.0 /im bands (Fig. [9l) are 
characteristic of low temperature, amorphous ices. Crys- 
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Figure 9. Profiles of the 3 /xm ice bands observed toward two 
background stars with the highest (solid and dotted lines) , and two 
background stars with the lowest (dashed and dash-dotted lines) 
extinctions and peak optical depths in the sample. The spectra 
are smoothed to a resolving power of 100. The profiles are very 
similar, indicating a lack of grain growth. 



talline ices would peak near 3.1 ^m. (jSmith et al.llT989l ). 

Finally, an absorptio n feature in t h e 5-8 /Ltm region 
referred to as C5 iii iBoogert et al.l (|2008D (see also 
iGibb fc Whittetl I2002D , has not been observed toward 
background stars (Table [5]). Although continuum reli- 
ability is generally a problem for this shallow feature, 
detections (> 3ct) were claimed in ^ 30% of the low 
mass YSOs and 50% of the massive YSOs. The fea- 
ture could be due to the formation of new molecules 
by energetic processing of the ices, possibly h ydrocar- 
bons ([Crecnberg et al. 1995) or ions (jSchutte fc Kham^ 
I2003f ). Its absence towards background stars is consistent 
with an origin in ices processed in the YSO environment, 
and not by external cosmic rays. 

5.3. Ices, Silicates, and the Extinction Curve 

Observations show that the properties of dust in dif- 
fuse and dense clouds are different. In contrast to diffuse 
clouds, the mid-infrared extinction in the dense cores 
studied here (^k >1 magn) stays hig h: at 25 ^,m. it 
is sti l l half of that a t 2 /im (Fig. ^3p; IChapman et al.l 
[2009t lMcClurel[2009l) . This is most likely caused by 
grain growth. Indeed, t he calculated i?v=5.5 ( "case B" ) 
extinction curve from [Weingartner fc Draind ()200lD is 
closer to the ob servations than their fly— 3. 1 "case A " 
model (Fig. H;; IChapman etld.l 120091: iMcChng [20091 . 
Nevertheless, even this extreme case of i?v=5-5 is sys- 
tematically too shallow at wavelengths above 5 /im, while 
the silicate features are too deep. The present work 
shows that th e contribution of ices, which are not in- 
cluded in the IWeingartner fc Draind (|2001f) models, to 
the extinction is not very large: ^25% in the 11-17 /J,m 
wavelength region. Instead, a model with larger grain 
sizes is needed. Such model would also need to take into 
account the observed weakness of the 9.7 /J,m band of 
sihcates (Figs.[2j: and [5^). 

Despite the need for larger grains in dense clouds, com- 
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pared to diffuse clouds, there is no evidence for progres- 
sive grain growth in dense clouds for Ak > 1 magn. 
The empirical extinction curve (Fig. [2]) satisfactorily fits 
most background stars in the sample, for the range of 
observed values of 1.3-6 mag^ The same conclusion 



can be drawn from Fig. 1 in iMcClurd ()2009f ). Also, 
the 3.0 /im band profiles are very similar over this A^ 
range (Fig. |9l), while a broadening and shift to longer 
wavelengths would be expected if grain growth occurred. 
Similarly, the observed linear relation of the 3.0 /im band 
with dust column density indicators such as the 9.7 /im 
band of silicates and Ak (Fig. [5|) indicates a lack of pro- 
gressive grain growth, since larger grains produce 3.0 /im 
bands with smaller peak optical depths. 

Recent work shows that the silicate band profiles 
are inconsistent with g rain growth beyond 1 //m 
(jvan Breemen et"al]|2011|). similar to w hat is found for 
the 3.0 jim ice band ([Smith et al1ll989f ). It is thus un- 
likely that grain growth can fully explain the shallower 
Ak versus rg.y relation in the dense medium compared 
to the diflFuse medium (Fig. [5^). Different dust composi- 
tions may play a role as well. The silicate bands in dense 
clouds show a short wavelength wing that could be due 
to a larger pyroxene/olivine abundance ratio ( §4.ip . For 
a more thorough invest igation of the Ak v e rsus T 9.7 re- 
lation one is referred to Ivan Breemen et all ()2011l ) . 

6. CONCLUSIONS AND FUTURE WORK 

The analysis of 2-25 /im photometry and spectra of 31 
stars tracing the quiescent medium in 16 isolated dense 
cores yields the following conclusions: 

1. An empirical "high resolution" extinction curve de- 
rived from one of the background stars confirms ex- 
tinction curves derived fro m broad band photome - 
try in high density clouds (jChapman et al.ll2009[ ). 
The curve remains remarkably flat above 10 /zm, 
such that at 25 /xm the extinction is still half that 
at 2.2 /xm. The main caveat is that it is assumed 
the near-infr ared cur v e is th e same as that derived 
bvUndcbcto uw et all ()2005[ ). 

2. The empirical extinction curve, in combination 
with model spectra of giants, and absorption spec- 
tra of H2O ice and silicates is able to fit the ob- 
served data of most background stars well. No sys- 
tematic deviations are seen as a function of core 
or Ak values (noting that the sample includes only 
high extinction lines of sight with Ak Jc^ 1-5 mag). 

3. Model fits of the attenuated stellar continuum are 
used to isolate the ice and dust absorption features. 
The 9.7 /zm band of silicates shows a shallower re- 
lation with Ak compared to the diffus e interstellar 
medium, thus confirming the results of iChiar et al.l 
(l2007h in many more lines of sight. Some back- 
ground stars, e.g., those behind the L328 core, fol- 
low the diffuse medium relation, however. 

4. The 3.0 /im band of II2O ice is detected in all ob- 
served sources. A strong correlation is observed 
between its peak optical depth and Ak- Its shape 
is independent of Ak, indicating a lack of grain 
growth at high Ak values. This may imply that the 



shallower relation of rg.y with Ak in the dense ISM 
compared to the diffuse ISM is not caused by grain 
growth. Possibly the dust composition is differ- 
ent, as also suggested by the broader 9.7 /im band 
in dense sight-lines (both YSOs and background 
stars). 

5. The 6.0 and 6.85 /im absorption bands are de- 
tected in most background stars. As for YSOs, 
the 6.0 /im band is only partially explained by the 
bending mode of H2O ice. Additional species, such 
as HCOOH, H2CO, CH3OH, and NH3 must be re- 
sponsible for this. Together with the consistent 
presence of the 6.85 /xm band, which may be due 
to NHJ, this shows that there is little difference be- 
tween the ices in quiescent core regions and those 
surrounding most YSOs. Most ices are apparently 
formed in the quiescent cloud phase. 

6. The discovery of solid CII3OH in the L-band spec- 
tra of several background stars strengthens the 

previous poin t. It confirms recent models by 

iCuppen et all (I2009D that CH3OH is efficiently 
formed on the grains in low temperature (T < 15 
K) clouds. Tight upper limits in previously studied 
sight-lines indicate that the ice composition varies 
in different environments at the onset of star forma- 
tion, however. For a proper comparison with the 
models, observed CII3OH/CO ratios are needed, 
i.e., more observations of solid CO toward back- 
ground stars. 

7. Signatures that are assigned to processing of the 
ices surrounding YSOs, such as profile changes in 
the 3 fim ice band, the shifted 6.85 /im band (i.e., 
a large C4/C3 ratio) and the presence of the broad 
C5 components, are not found toward the back- 
ground stars. 

The present work studies ices and the extinction 
curve in the quiescent medium in isolated dense 
cores. YSOs in these cores are not well studied. 
On the other hand, the ices are well studied to- 
ward YS Os in large clouds (Taurus, Serpens, Oph 
Perseus; [Boogert et al.l 120081: iPontoppidan et all 120081: 
lOberg et all 120081 : iBottinelh et~ani2010D. but th e back- 
ground star samples are small (jKnez et al.ll2005l ). Thus, 
future studies will need to compare YSOs and back- 
ground stars in the same environments to address the 
origin of the observed variations of the CII3OH abun- 
dance and the strength of the 6.0 and 6.8 fim ice bands 
(Figs. |3]and[7]). Furthermore, while there is now much 
observational evidence that the dust properties in the 
dense and diffuse medium are quite different, the cause 
of this is not understood, warranting further theoretical 
and perhaps laboratory investigations. 
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